GSMN-ML- a genome scale metabolic network reconstruction of the obligate human pathogen Mycobacterium leprae.

Leprosy, caused by Mycobacterium leprae, has plagued humanity for thousands of years and continues to cause morbidity, disability and stigmatization in two to three million people today. Although effective treatment is available, the disease incidence has remained approximately constant for decades so new approaches, such as vaccine or new drugs, are urgently needed for control. Research is however hampered by the pathogen's obligate intracellular lifestyle and the fact that it has never been grown in vitro. Consequently, despite the availability of its complete genome sequence, fundamental questions regarding the biology of the pathogen, such as its metabolism, remain largely unexplored. In order to explore the metabolism of the leprosy bacillus with a long-term aim of developing a medium to grow the pathogen in vitro, we reconstructed an in silico genome scale metabolic model of the bacillus, GSMN-ML. The model was used to explore the growth and biomass production capabilities of the pathogen with a range of nutrient sources, such as amino acids, glucose, glycerol and metabolic intermediates. We also used the model to analyze RNA-seq data from M. leprae grown in mouse foot pads, and performed Differential Producibility Analysis to identify metabolic pathways that appear to be active during intracellular growth of the pathogen, which included pathways for central carbon metabolism, co-factor, lipids, amino acids, nucleotides and cell wall synthesis. The GSMN-ML model is thereby a useful in silico tool that can be used to explore the metabolism of the leprosy bacillus, analyze functional genomic experimental data, generate predictions of nutrients required for growth of the bacillus in vitro and identify novel drug targets.

Network analysis has diverse roles in drug discovery.

Computational biologists use network analysis to uncover relationships between various data types of interest for drug discovery. For example, signalling and metabolic pathways are commonly used to understand disease states and drug mechanisms. However, several other flavours of network analysis techniques are also applicable in a drug discovery context. Recent advances include networks that encompass relationships between diseases, molecular mechanisms and gene targets. Even social networks that mirror interactions within the scientific community are helping to foster collaborations and novel research. We review how these different types of network analysis approaches facilitate drug discovery and their associated challenges.

¹³C metabolic flux analysis identifies an unusual route for pyruvate dissimilation in mycobacteria which requires isocitrate lyase and carbon dioxide fixation.

Mycobacterium tuberculosis requires the enzyme isocitrate lyase (ICL) for growth and virulence in vivo. The demonstration that M. tuberculosis also requires ICL for survival during nutrient starvation and has a role during steady state growth in a glycerol limited chemostat indicates a function for this enzyme which extends beyond fat metabolism. As isocitrate lyase is a potential drug target elucidating the role of this enzyme is of importance; however, the role of isocitrate lyase has never been investigated at the level of in vivo fluxes. Here we show that deletion of one of the two icl genes impairs the replication of Mycobacterium bovis BCG at slow growth rate in a carbon limited chemostat. In order to further understand the role of isocitrate lyase in the central metabolism of mycobacteria the effect of growth rate on the in vivo fluxes was studied for the first time using ¹³C-metabolic flux analysis (MFA). Tracer experiments were performed with steady state chemostat cultures of BCG or M. tuberculosis supplied with ¹³C labeled glycerol or sodium bicarbonate. Through measurements of the ¹³C isotopomer labeling patterns in protein-derived amino acids and enzymatic activity assays we have identified the activity of a novel pathway for pyruvate dissimilation. We named this the GAS pathway because it utilizes the Glyoxylate shunt and Anapleurotic reactions for oxidation of pyruvate, and Succinyl CoA synthetase for the generation of succinyl CoA combined with a very low flux through the succinate--oxaloacetate segment of the tricarboxylic acid cycle. We confirm that M. tuberculosis can fix carbon from CO2 into biomass. As the human host is abundant in CO2 this finding requires further investigation in vivo as CO2 fixation may provide a point of vulnerability that could be targeted with novel drugs. This study also provides a platform for further studies into the metabolism of M. tuberculosis using ¹³C-MFA.

Differential producibility analysis (DPA) of transcriptomic data with metabolic networks: deconstructing the metabolic response of M. tuberculosis.

A general paucity of knowledge about the metabolic state of Mycobacterium tuberculosis within the host environment is a major factor impeding development of novel drugs against tuberculosis. Current experimental methods do not allow direct determination of the global metabolic state of a bacterial pathogen in vivo, but the transcriptional activity of all encoded genes has been investigated in numerous microarray studies. We describe a novel algorithm, Differential Producibility Analysis (DPA) that uses a metabolic network to extract metabolic signals from transcriptome data. The method utilizes Flux Balance Analysis (FBA) to identify the set of genes that affect the ability to produce each metabolite in the network. Subsequently, Rank Product Analysis is used to identify those metabolites predicted to be most affected by a transcriptional signal. We first apply DPA to investigate the metabolic response of E. coli to both anaerobic growth and inactivation of the FNR global regulator. DPA successfully extracts metabolic signals that correspond to experimental data and provides novel metabolic insights. We next apply DPA to investigate the metabolic response of M. tuberculosis to the macrophage environment, human sputum and a range of in vitro environmental perturbations. The analysis revealed a previously unrecognized feature of the response of M. tuberculosis to the macrophage environment: a down-regulation of genes influencing metabolites in central metabolism and concomitant up-regulation of genes that influence synthesis of cell wall components and virulence factors. DPA suggests that a significant feature of the response of the tubercle bacillus to the intracellular environment is a channeling of resources towards remodeling of its cell envelope, possibly in preparation for attack by host defenses. DPA may be used to unravel the mechanisms of virulence and persistence of M. tuberculosis and other pathogens and may have general application for extracting metabolic signals from other "-omics" data.

The application of next-generation sequencing technologies to drug discovery and development.

Next-generation sequencing (NGS) technologies represent a paradigm shift in sequencing capability. The technology has already been extensively applied to biological research, resulting in significant and remarkable insights into the molecular biology of cells. In this review, we focus on current and potential applications of the technology as applied to the drug discovery and development process. Early applications have focused on the oncology and infectious disease therapeutic areas, with emerging use in biopharmaceutical development and vaccine production in evidence. Although this technology has great potential, significant challenges remain, particularly around the storage, transfer and analysis of the substantial data sets generated.

The genetic requirements for fast and slow growth in mycobacteria.

Mycobacterium tuberculosis infects a third of the world's population. Primary tuberculosis involving active fast bacterial replication is often followed by asymptomatic latent tuberculosis, which is characterised by slow or non-replicating bacteria. Reactivation of the latent infection involving a switch back to active bacterial replication can lead to post-primary transmissible tuberculosis. Mycobacterial mechanisms involved in slow growth or switching growth rate provide rational targets for the development of new drugs against persistent mycobacterial infection. Using chemostat culture to control growth rate, we screened a transposon mutant library by Transposon site hybridization (TraSH) selection to define the genetic requirements for slow and fast growth of Mycobacterium bovis (BCG) and for the requirements of switching growth rate. We identified 84 genes that are exclusively required for slow growth (69 hours doubling time) and 256 genes required for switching from slow to fast growth. To validate these findings we performed experiments using individual M. tuberculosis and M. bovis BCG knock out mutants. We have demonstrated that growth rate control is a carefully orchestrated process which requires a distinct set of genes encoding several virulence determinants, gene regulators, and metabolic enzymes. The mce1 locus appears to be a component of the switch to slow growth rate, which is consistent with the proposed role in virulence of M. tuberculosis. These results suggest novel perspectives for unravelling the mechanisms involved in the switch between acute and persistent TB infections and provide a means to study aspects of this important phenomenon in vitro.

Selection of objective function in genome scale flux balance analysis for process feed development in antibiotic production.

Using flux variability analysis of a genome scale metabolic network of Streptomyces coelicolor, a series of reactions were identified, from disparate pathways that could be combined into an actinorhodin-generating mini-network. Candidate process feed nutrients that might be expected to influence this network were used in process simulations and in silico predictions compared to experimental findings. Ranking potential process feeds by flux balance analysis optimisation, using either growth or antibiotic production as objective function, did not correlate with experimental actinorhodin yields in fed processes. However, the effect of the feeds on glucose assimilation rate (using glucose uptake as objective function) ranked them in the same order as in vivo antibiotic production efficiency, consistent with results of a robustness analysis of the effect of glucose assimilation on actinorhodin production.

GSMN-TB: a web-based genome-scale network model of Mycobacterium tuberculosis metabolism.

BACKGROUND: An impediment to the rational development of novel drugs against tuberculosis (TB) is a general paucity of knowledge concerning the metabolism of Mycobacterium tuberculosis, particularly during infection. Constraint-based modeling provides a novel approach to investigating microbial metabolism but has not yet been applied to genome-scale modeling of M. tuberculosis. RESULTS: GSMN-TB, a genome-scale metabolic model of M. tuberculosis, was constructed, consisting of 849 unique reactions and 739 metabolites, and involving 726 genes. The model was calibrated by growing Mycobacterium bovis bacille Calmette Guérin in continuous culture and steady-state growth parameters were measured. Flux balance analysis was used to calculate substrate consumption rates, which were shown to correspond closely to experimentally determined values. Predictions of gene essentiality were also made by flux balance analysis simulation and were compared with global mutagenesis data for M. tuberculosis grown in vitro. A prediction accuracy of 78% was achieved. Known drug targets were predicted to be essential by the model. The model demonstrated a potential role for the enzyme isocitrate lyase during the slow growth of mycobacteria, and this hypothesis was experimentally verified. An interactive web-based version of the model is available. CONCLUSION: The GSMN-TB model successfully simulated many of the growth properties of M. tuberculosis. The model provides a means to examine the metabolic flexibility of bacteria and predict the phenotype of mutants, and it highlights previously unexplored features of M. tuberculosis metabolism.